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Abstract. We show existence and uniqueness for the solutions of the regularity and the 

Neumann problems for harmonic functions on Lipschitz domains with data in the Hardy 

2 

spaces H p (dD)(H p (dD)), p > — — e, where D CM 2 and e is a (small) number depending on 

o 

the Lipschitz nature of D. This in turn implies that solutions to the Dirichlet problem with 
data in the Holder class C^ 2+£ (dD) are themselves in C 1 / 2+£ (I?). Both of these results 
are sharp. In fact, we prove a more general statement regarding the H p solvability for 
divergence form elliptic equations with bounded measurable coefficients. 

We also provide H 2 ^ 3 ~ E and C x / 2+e solvability result for the regularity and Dirichlet 
problem for the biharmonic equation on Lipschitz domains. 



1. Introduction and Main Results 

In this paper we study the Dirichlet and Neumann problems for harmonic functions on 
Lipschitz domains and their biharmonic counterparts. More precisely let X, Y, Z be function 
spaces on the boundary dD of D. Then 



(D 



Xj 



Au 

u\dD 
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= 
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is the Dirichlet problem with underlying space X, and 

An =0 

U1S dD 

e Y 



dD 

M(Vu) 



Au 

u\dD 

M(Vu) 



= 
= h 
e Z 



du 

are the Neumann and regularity problems. Here iV is the outer normal vector to D, — — = 

(N, Vw) and M(u) is the usual non-tangential maximal function of u. 

In this setting the canonical choices are X = Y = L p (dD), Z = L\{dD)- the space of 
functions with one distributional derivative in L p (dD). In the sequel, we will slightly abuse 
notations by using D p instead of -D_z>, N p instead of Nlp etc. 

1.1. Harmonic functions. We state now the classical results related to the IP theory. 



Theorem [3, 12, 4, 5] Let D C 



be a connected Lipschitz domain. Then 
l 
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1. there exists an e = e(D) > such that for 2 — e < p < oo and / G L p (dD) there is an 
unique solution to _D P . Moreover, there is the apriori estimate 
\\M(u)\\ LP{dD) <C\\f\\ LP{dD) . 

2. there exists ane = s(-D) > such that for 1 < p < 2 + e and g G L p (dD), j QD gda = 
there is an unique (up to a constant) solution to N p . There is the apriori estimate 
\\M(Vu)\\ LP{dD) <C\\g\\ LP{dD) . 

3. there exists ane = e(D) > such that for 1 < p < 2 + e and h G L\{dD) there is an 
unique solution to R p . There is the apriori estimate 
\\M(Vu)\\ LP{dD) <C\\h\\ LP{dD) . 

This theorem summarizes the results in [3], but some earlier version and ideas originated 
in [12]. Actually,the LP theory described above is a consequence of the duality between the 
Dirichlet and regularuty problems, the L 2 solvability for all three problems and the following 
endpoint result due to Dahlberg and Kenig. 



Theorem 1 (Dahlberg-Kenig) 

Then 

1. 



Let D C W 1 be a connected star-like Lipschitz domain. 



Au 
du 

M(Vu) 



= 

= / 
GL 1 



is uniquely solvable provided f G H 1 (dD) and \\M(Vu)\\ Ll ^ dD ^ < C\\f\\ Hl ^ dD y 
2. Given f G H\{dD) there exists an unique solution to 



Au 

u\dD 

M(Vu) 



= 

= / 
GL 1 



Moreover \\M(Vu)\\ L1(dD) < C\\f\\ Hl(dD) . 



As a corollary to this result, one proves (weak) maximum principle, solvability for BMO 
data etc. We refer to [3] for excellent treatise of these questions. 

Recently, Brown [1] was able to extend Theorem 1 to show that there exists e = e(D) , 
such that for 1 — e < p < 1 the Neumann problem N p is still uniquely solvable with the usual 
estimates ||M(Vu)|| < C||/||# P - This result has the interesting corollary that the double 
layer potential is invertible operator on the Holder space C a (dD) for a close to zero and 
thus we have a representation formula for the solutions of the Dirichlet problem with C a 
data. This raises the following natural question, see Question 3.2.10 in [5]. 

Question 1. Are R p and N p solvable for p significantly below one? What does that imply 
for solutions of the Dirichlet problem with C a data for a significantly above zero? 

The purpose of this paper is to establish the optimal p range for solvability of both Rp 
and N p in dimension two. That is our Theorem 2 below. Let us remark only that known 
counterexamples in dimensions bigger than two imply that the Neumann problem N p may 
not be solvable for p < 1 — e(D), i.e. for fixed p < 1, there exists a Lipschitz domain D such 
that N p (D) is not uniquely solvable. 
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Theorem 2. Let D C M, 2 be a star-like Lipschitz domain with connected boundary. There 
exists e = e(D), such that for 2/3 — e < p < 1 and < a < 1/2 + e 

1. The Neumann problem 

Au = 
O^ldD - J e H [c)V) 

is uniquely solvable and \\M(Vu)\\ LP ^ dD ^ < C\\f\\ HP ^ 9D y 

2. The regularity problem 

Au = 



u\ 9 d = fe H{{dD) 



has unique solution and \\M(Vu)\\ LP ^ dD ^ < C\\f\\ H p^ dD y 
3. The Dirichlet problem 



(D a ) 



Au = 
u\ 9D = fe C a (dD) 

has unique solution and \\u\\ Ca ^ < C\\f\\ Ca , dD y 
Moreover, the ranges 2/3 — e < p and a < 1/2 + e are sharp. 

In fact, we consider more general divergence form elliptic equations in the form 
div(A(Vu)) = 0, where A is a symmetric, elliptic matrix with real-valued bounded measur- 
able coefficients. We prove the following theorem. 

/ a (x) b(x) \ 

Theorem 3. Let A(x,t) = A(x) = ( ^) I )( J be a real, symmetric, uniformly elliptic 

matrix with bounded and measurable coefficients, independent of the time variable. Then 
there exists e = e(D) > such that for 2/3 — e < p < 2 + e 

1. the Neumann problem in the upper half-space = {(x, t) : t > 0} 

div(AVw) = for t > 0, 
AVu(x, 0) • (0, —1)| = f E H P (M}) 

has unique solution and ||M(Vu)|| iP < C'H/H^p. 

div(AVu) = for t > 0, 
«(a;,0) = he Hf{R r ) 

has unique solution and \\M(Wu)\\ LP < C\\h\\ H p. 

Moreover, the range 2/3 — e < p is sharp. 

Remarks 

• For divergence form equations div(AVw) = 0, A = A(x, t) one cannot expect solvability 
even for D 2 or N 2 . Indeed, counterexamples show that unless we require radial inde- 
pendence for such a problem in the unit ball, we may encounter non-uniqueness for N 2 , 
see [6] and [5], p. 63. 

We will however work in the upper half-space instead of the unit ball. These two prob- 
lems are not so much different. In fact, the appropriate assumption in the upper-half 
space is time independence (see [5], p. 68 for a relevant discussion) and so our theorem 3 



(N P ) 



2. The regularity problem in the upper half-space m " 2 

(R P ) 
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is formulated in that fashion. Let us remark only, that the problems D 2 , R 2 and N 2 are 
all solvable for matrices A = A(x) with time independent coefficients ([6]). We make 
use of these facts later on in our proofs. 
• The restriction to the upper-half space in theorem 3 is just for technical reasons. In 
fact, one can state the theorem for a general Lipschitz domain in M 2 . The following 
argument shows that for the Dirichlet (regularity) problem. 

Let D be the domain above the Lipschitz graph t = tp(x) and u solves the Dirichlet problem 
div(A(x)Vu) = 0, u(x,tp(x)) = f(x). Define t) = (x, t — ip(x)) and set u($(x,t)) = 
u(x,t). It is not difficult to check that u : — > M 1 is a solution to div(AVu) = 0, 
u(x, 0) = /, where 



A(x) 



A(x) 



V(x) 
1 



1 

-ip'(x) 1 

In particular, we have shown that Theorem 3 implies parts one and two of Theorem 2. 

Unfortunately, at this moment we cannot claim part three of our Theorem 2 for general 
divergence form elliptic equations with time independent coefficients. Our proof for harmonic 
functions is based on Brown's duality technique for the double-layer potential, which does 
not seem to generalize in the setting of Theorem 3. Thus we pose the following: 

Question 2. Assume that A = A(x) is a real, symmetric elliptic matrix. Prove that the 
Dirichlet problem in the upper half-space M. 2 , 



div(AVu) = for t > 0, 
u\ dD = fe C a (R r ) 

for a < l/2 + e(A). 



(D a ) 

is solvable with \\u\\ Ca ^ R ^ < C\\f\\ c , 

We now state the following result, which gives a connection between the Neumann Hardy 
spaces and the usual atomic Hardy spaces. This is an extension of Theorem 2.3.18 in [5], for 
the case p < 1. 

Theorem 4. Let D c8 2 is a domain above Lipschitz graph and u is a harmonic function 
on D that satisfies u\qd G Hf(dD), 2/3 — e < p < 1. Then 

du 

'\ u \\m(dDy 



dN 



< 



HP(dD) 



du 



Conversely, given f G H p (dD), there exists a harmonic function u, such that = f and 



HP(dD) ' 



1.2. Biharmonic functions. For the biharmonic equation, we consider the Dirichlet prob- 
lem 

A 2 u 



BD r 



U\dD 

du 

\dD 



= 

= /o 



= E fjNj, 
3=1 
< 00, 



dN 

\\M(Vu)\\ LP{dD) 

where Ni, N 2 , . . . ,N n are the components of the normal vector and /o, /1, /2, • • ■ , fn satisfy 
the compatibility condition (/ , /1, / 2 , . . . , /„) G WA 2 (dD) (cf. [10]). The regularity problem 



BR P 
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is 

A 2 u =0 
D 2 u\ dD = f 

3=1 

\\M(VVu)\\ LP{9D) <oo 

The L 2 theory for these problems (with the necessary adjustments for the order of the 
derivatives) is very similar to the harmonic case and we present it in Section 2. We have the 
following results in two dimensions. 

Theorem 5. There exists an e = e(D) > 0, such that if < a < 1/2 + e, fi, f 2 G 
C a (dD) fl L 2 (dD), then the unique L 2 solution to BD a satisfies Vu G C a (D). In fact, 

2 

\\Vu\\ Ca{D) + sup dist(X, dD)- l - a \u{X) - u(X*) -(X- X*, Vu(X*))\ < C V \\fj\\ Ca , D) , 

- ' xeD v ; 

3=1 

where C is a constant depending only on the Lipschitz nature of D and X* is the projection 
of X along the "time" axis onto dD. Moreover the range a < 1/2 + e is sharp. 

Theorem 6. There exists an e = s(D) > 0, such that the regularity problem BR p with 
2/3 — e < p < 2 + e, (/, g) G Hf(dD) x H p (dD) has unique (up to a constant) solution. 
Moreover the estimate 

\\M(WVu)\\ LP{dD) < c\\v T J\\ HP{9D) + \\g\\ HP{dD) , 

holds with a constant C depending only on the Lipschitz nature of D. The range 2/3 — e < p 
is sharp. 



2. Preliminaries 

We separate this section into two parts - about harmonic and biharmonic functions re- 
spectively. The corresponding equations exhibit some common features like the L 2 theory, 
but there are some dissimilarities as well. We try to present the similarities in the technically 
simpler harmonic context and we briefly outline some specifics for the biharmonic operator. 
We constantly refer in the text to the papers [9], [10] for the necessary background results. 
Since the two dimensional case is of utmost interest to us, we sometimes avoid the explicit 
formulas (with the inevitable technicalities that arise) for dimensions higher than two. 



2.1. Harmonic functions. Let D C W 1 be a Lipschitz domain, such that D,D C are con- 
nected. For technical reasons, we restrict our attention to the case of domains above Lipschitz 
graphs, i.e. 

D = {{x, t) : t > (p(x)}, ip : M n_1 -> K 1 , 
\(p(x) - <p(y)\ < M\x - y\. 



The surface measure on dD is defined via the usual da = a/1 + \V(f\ 2 dx. 
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Following [3], we introduce the atomic Hardy spaces H p (dD) for 1 > p > (n — l)/n. First 
an H p (dD) atom is a function a : dD — > R, such that 

supp(a) C B(Q, d) = {P G dD : \P - Q\ < d}, 
J ada = 0, 

Hall , < Cd (n - 1)(1/2 ~ 1/p) 

\\ a \\L 2 {dD) — ° tt 

Then, 

H p (dD) = {Y J \a t :Y J \M P < °c} 



i/p 



where aj are H p (dD) atoms. 

We also define H p (dD) atoms by requiring that 

supp(a) C B(Q, d) = {PedD:\P-Q\<d}, 

(d dtp d \ 
— h — — — — J -u. The space Hf(dD) of distributions with one 

derivative in H p (dD) may be defined as the V span of such atoms. It is well known that 
H 1 (dD) C L l (dD), while the spaces H p (dD), p < 1 contain non-integrable distributions. 
Sometimes, we will abuse notations by writing a(x), instead of a(x,ip(x)). Observe that 
d 

V Ti a = —a(x,ip(x)). 

We also define the (homogeneous) Holder spaces C a (dD), < a < 1 by 

C a (0D) = j/ : dD - M 1 : ||/|| Ca(aD) = sup ~ } . 

Ti — 1 

We remark that the Holder spaces C a and the Hardy spaces H p , p = can be 

n — 1 + a 

paired in the sense that every element in one of them defines via integration a continuous 
linear functional on the other (cf. [11] p. 130). Let 

|^|2— n 



T(x) = I (™-2R 



n > 2 



— In x n = 2 
2tt 1 1 

be the fundamental solution for the Laplace's equation in R™. Define the single and double 
layer potentials S and /C by 

S(f)(X) = p.v. Jr(X-Q)f(Q)da(Q) : xeR n \dD 

dD 

}C(f)(X) = p.v. J J?jL(X-Q)f(Q)da(Q), xeW\dD 



3D 
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We also define the formal adjoint of /C 



K*(f)(X) = p.v. J - Q)f(Q)da(Q). 



dD 

2.2. Biharmonic functions. We start with the L 2 theory for the biharmonic equation, due 
to Kenig and Verchota [7] (see also Theorem 3.7 in [9]). 

Proposition 1. Let D dW 1 be a Lipschitz domain. Then there exists an e > depending 
only on the Lipschitz nature of D, such that forp : 2 — e < p < 2 + e the equation 

A 2 u = 

u\dD = f 

{N,Vu)\ dD =g 
\\M(Vu)\\ LP{dD) <oo 

is uniquely solvable. In addition, there are the estimates 

\Vu(X)\ < dist(X,9 J D)-("- 1 )/f, 
\\M(Vu)\\ LP9D) <\\Vu\\ dD . 

There is a also the regularity result, which we now state( cf. Theorem 4.6, [9]). 

Proposition 2. Let D C lR n be a Lipschitz domain. Then there exists an e > depending 
only on the Lipschitz nature of D, such that for p : 2 — e < p < 2 + e the equation 

A 2 u =0 

D n u\dD = f 

ZU(V Tj ,VD 3 u)\ dD =g 
||M(VVu)|| 

LP (3D) < 00 

is uniquely solvable. In addition, there are the estimates 
• |VVu(X)| < dist(X,9 J D)-( n - 1 )/ p , 



ll^vv^ii^^E.dlVT./ll 



LP(dD) 



+ \\g\ 



LP {3D) 



The fundamental solution of the biharmonic equation in two dimensions is 

Based on the L 2 theory, one is able to define the Green's function as follows. Let X be inside 
the domain and fix a point X dD. Let 

E(X -Q)- S(X - Q) = ^-(\X - Q\ 2 In \X - Q\ - \X - Q\ 2 In \X - Q\), 



8tt 



HQ) = 

f 3 (Q) = Dj(f (Q)), je{l,2}. 
Consider then the unique solution 7x00 to the problem 

A 2 u = 0, 

u\dD = fo, 

Vm| 9 D ={/l,/2}- 



\L 2 (dD) 



< 1. 



Since V^/j ~ \X — Q\ we have by the L 2 regularity result that ||M(VV7x) 
The Green's function can be defined as G(X, Y) = E(X -Y)- S(X -Y)- j x (Y). Observe 
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that since two tangential derivatives of the explicit quantity S(X — Y) — S(X — Y) also 
belong to L 2 (dD), we can conclude 

(1) \\M(VVG(X,.))\\ L2{dD) <l. 

Later on, we will be able to show a much stronger estimate than (1) when the integration is 
over dyadic pieces away from the origin. 

Integration by parts shows that one has the following representation formula for bihar- 
monic functions (cf. [10]): 

u(Q)—A Q G(X,Q) + / —-(Q)AqG(X, Q). 

dD dD 

3. Existence and Uniqueness for harmonic functions 

In this section, we prove existence and uniqueness for R p and N p , and Dirichlet problem 
with Holder data in Theorem 3. Our plan is as follows. First, we show that the regularity 
and Neumann problems are equivalent, i.e. if one can solve the regularity problem uniquely 
in Hf(dD), then one can solve uniquely the Neumann problem in H p (dD) and vice versa. 
Secondly, we show that the R p is solvable in Hf(dD), p > 2/3 — e and finally we prove 
uniqueness for R p . The uniqueness result will be almost automatic in view of Lemma 2.2 in 
[1] and the usual L q uniqueness result for the Dirichlet problem, q > 2 — e. Finally, we show 
the existence and uniqueness result for the Dirichlet problem with C a (dD) data in Theorem 
2. 

3.1. Equivalence for the regularity and Neumann problem. Let u satisfy the Neu- 
mann problem N p with data / e H p (dD). By the properties of H p (dD), we consider without 

t 

loss of generality only smooth compactly supported data /. Let U-i(x, t) — J u(x, z)dz. De- 

o 

fine 

v = a(x)(u-i) x + b(x)(u-i) t . 

It is not difficult to check that 



v x (x, 0) = (a(u_i)a. + b(u-i) t ) x = -(&(u-i)irt + c(u_i) tt ) = -bu x - cu t = f(x). 
We observe that v satisfies 

' div(iVw) = 
v(x,0) =g(x)' 



where A = — and g is the (unique) function with g'ix) — fix), f g — 0. 

ac — b 2 

Since A is an uniformly elliptic matrix with time independent coefficients, so is A. Note 
that starting with a solution of a Neumann problem, we have produced an L 2 solution 

to an associated regularity problem. Moreover, since Vi> =1 , J Vm, one obtains 

pointwise equivalence M(Vv ) ~ M(Vu). Hence, if one can prove estimates for R p 

\\M(Vv)\\ LP{dD) <\\g\\ Hf{dD) , 
they would imply the corresponding estimates for N p 

\\M(Vu)\\ LndD) <\\f\\ HP(9D) . 
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Thus, we have showed that if one can solve the regularity problem in Hf(dD), then one can 
also solve the Neumann problem. The reverse implication can be proved by retracing back 
the argument above, so we omit the details. 

We note that the equivalence of the regularity and Neumann problems in the sense de- 
scribed above is purely two dimensional phenomena. Actually, in the important case of 
the Laplace's equation, it is not difficult to check that u and v above are in fact conjugate 
harmonic functions and thus one cannot expect the equivalence to persist in higher dimen- 
sions. Actually, by the equivalence of the regularity and Neumann problem and the existence 
results of Theorem 2 (to be proved below), we establish Theorem 4. 

3.2. Solvability for the regularity problem in Hf(dD). By well known approximation 
techniques (see for example [5], section 1.10), it will suffice to prove the estimate 

(3) \\M(Vu)\\ LP <C\\g\\ HPi 

for solutions u of R p corresponding to smooth matrices A and smooth data g, which are 
known to exist, as long as the constant C is independent of everything, but the ellipticity 
constant of A. 

Thanks to the "atomic" nature of iJ^IR 1 ), one can take g to be a iJf-atom. Indeed, if we 
assume (3) for atoms and take into account the p-subaditivity of the LP quasi- norm (p < 1) 
we get for g = £ 

\\M(vu g )\\ p p < J2 \mM(vu gi )\\ p p < E m p s imiv 

Simple dilation and translation argument allows us to reduce to the case of an unit atom, 
i.e. 

1. supp g C [-1, 1], 

2. IbIL, ML £ !■ 

For r e (1/2, 1), consider the intervals R] = (2 j r, 2 j+1 /r) and Rj := R). Let 

( 100 (2* t -x) x< 2-V, 

gj(x) = < 2 j r<x< 2^' +1 /r, 

[ 100(x-2^ +1 /r), x>T +1 /t, 

and Vt] = {(x,t) : t > q](x)}. Ob serve that since QZ is a domain above Lipschitz graph, the 
L 2 theory for divergence form equations with time independent coefficients applies to it (see 
the discussion after Theorem 3). We have the following lemma. 

Lemma 1. Let u be the unique L 2 solution to the problem 

div(AVw) = 
u(x,0) =g(x), 

where g is an unit atom in LLfiW 1 ). Then there exists e > 0, such that 



J M(Vu) 2 < C £ 2^ £ - 2 ^. 
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Let us show that Lemma 1 implies (3). By Holder's inequality, 

/ \ P/2 



J M{Vuf < \R J \ 1 - p/2 J M{Vu) 

Rj \Rj J 



< 2 J'(l-p/2+p(-e-2)/2) 



and for every p > 2/3 — 0(e), the series J2j Jr M(Vu) p converges. 
Thus, it remains to prove Lemma 1. 

Proof. (Lemma 1) We use the standard Cacciopoli type argument. By the L 2 theory for f2J 
and since Rj C <9f2J, we derive 

i 



iVul 



(4) Jm(Vu) 2 <J J M(Vu) 2 dr<2-> J 

R 3 1/2 dttj 0, 

1/2 

where flj = Vt- . Break Qj into "good" and "bad" part, so that 

Gj = Qj n{t> 2 j } 

oo 

On the good part, we further decompose Gj = \J Gj, so that Gj = Gj n {t ~ 2 j+k }. For 

k=i 

each Gj, one applies the usual interior estimates for the solution (cf. [5]), to get 

/ ,„ ,9 , ^— <>(i-i-k\ I , 

\U\ 



(5) J \Vu\ 2 < 2- 2 ^ J 



By the D 2 solvability, we can always estimate |M| L 2( G fc) < 2^ +k ^ 2 , which would give the 
desired estimate, except fot the extra decay factor 2~ £j . 

For the "bad" part, select an even function i/j e C^(R 2 ), so that supp ip C (1/4,4) x (0,2) 
and ip{x,t) = 1 for G (1/2,2) x (0,1). Denote ipj := ip(2~ j -). Observe that since 

u(x, 0) = on Rj, we may extend u(x, t) for t < across Rj as an even function. By the 
ellipticity of A and the divergence theorem, one then derives 

(6) J \Vu\ 2 < J(AV(u^),V(u^))<2-> J |V«||«|^-<2- J '||V«|| L2(c . ) ||«|| La(c . ) , 

Bj R2 R 2 

where Cj D Bj is again a box with sides ~ 2K 

From (6), R 2 , D 2 solvability and by iterating (6) (we will get back to this point later on), 
we easily get the bound f n |Vw| 2 < Cs2~ j 2 Sj for all positive 5. This estimate, together with 

(4) imply Lemma 1 without the crucial term 2~ £J . 

The usual approach to get the improvement 2~ £j is to use Sobolev embedding if 1 (IR n ) C 
L 2n ^ n ~ 2 ^ (R n ) , which unfortunately fails for dimension two. We use instead the following 
multiplicative variant of Sobolev embedding 

(7) IMIl4 ( r2) < l|M|li2 2 (R2) ||VM||^ 2 
We have the following proposition. 
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Proposition 3. Suppose 

(8) HVu||£ a( Q.) <1- 3 \\u\\ L2{Pj) \\Vu\\ L2{P . v 

(9) \\M(u)\\ L2 . S{Rl) <l, 

(10) ||M(V«)|| L2(K1) <1. 

where Qj C Pj are boxes with sides ~ 2 J . Then there exists e = 0(5) > 0, such that 

(11) !lv«lb (Qi) < 2(- 1 / 2 -^'. 

It is clear that a direct application of (11) gives the estimate for the "bad" part, while 
for the good part one applies (11) for G'j and summation in k > then gives (5). Hence, to 
complete the proof of Lemma 1, it remains to prove Proposition 3. □ 

Proof. (Proposition 3) 

Apply (7) for u(x,t)ipj(x — xo,t — to)) where (xo,to) are suitably chosen so that ipj(x — 
xo, t — t ) = 1 on Pj and supp ipj(- — xq, ■ — t ) C 4Pj. Cauchy-Schwartz and (9) yield 



/ 



u? < \\u\\ l2 -s {Pj) \\u\\ l(2 . sY{Pj) < 2»VM-^) 2 W-*)'|| u || L4(pj) < by (7) 



< (1/2+1/(2-5)' ) / 1 1 1 1 1/2 ||v ,,1/2 +0-i/H v \\ \< 

< 2 ,(l/2 + l/(2- < 5)') (2J /4 ||Vw|| l/2 4 ^ ) + 1} < 2 ( 5/4 -0(^ ||W|| l/2 4 ^ ) + 2j2 - 0(5)j - 

where 0(5) is a positive number of the order of 5. From the preceding estimate and (8), we 
get 

(12) / |V M | 2 < 2-^°^u\\%\ 4Pj) + 2-/-^||V.|| L2(4Pj) 

We can perform now the following iteration procedure. Call \j = || Vm|| l2 (q ) and fij = 
\\Vu\\ L 2^ p .y Clearly (12) reads as 

(13) \) < 2-^2-°^ /if + 2-^ 2 2-°W> J , 

Since (10) allows us to bound fij < C2^ 2 , one gets from (13) improvement for the bounds 
for \j, jij. We continue in that fashion and use the improved bounds back at (13). That way 
one gets an improvement at every step. One has 

for some e = 0(5) and the proof is complete. □ 

3.3. Uniqueness for the regularity problem in if*. The uniqueness result is almost 
automatic in two dimension due to the following lemma of Brown [1] , which we state verba- 
tim. 

Lemma 2. Let D C W 1 be a connected Lipschitz domain and suppose that u is harmonic 
in D. Let X* be a fixed point in D and suppose u(X*) = 0. For p < n — 1 and p* = 
(n — l)p/ (n — 1 — p) we have 

(14) \\M(u)\\ LP * {9D) <C\\M(Vu)\\ LP{9D) , 
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where the constant C depends only on the distance of X* to the boundary, p and the Lipschitz 
character of dD. 

Carefull inspection of the proof shows that one can relax the harmonicity assumptions 
on u, by requiring that u satisfies a divergence form equation. Indeed in the proof of (14), 
Brown uses interior estimates and the equivalence of the square function with the non- 
tangential maximal function in L 2 , which are of course available for solutions of divergence 
form equations as well. 

In contrast with the higher dimensionional case, where one needs to have an additional 
argument to prove uniqueness for N p , p > 1 — e (cf. [1]), the two dimensional uniqueness 
result follows from the Brown's lemma for the range 1 > p > 2/3 — e and uniqueness for D q , 
q > 2—e. To this end, assume that u solves R p with zero data, such that ||M(Vu) || iP ( K i) < oo. 
From (14) we get 

\\M(u)\\ LP * {Rl) <C\\M(Vu)\\ LP{Rly 

Observe that since 1 > p > 2/3 — e, we have 2 — 0(e) < p* = p/(l — p) < oo and therefore 
u solves a Dirichlet problem (with zero data), with ||M(u)|| LP ^ E i^ < oo. Thus u = by the 
uniqueness for D p *. 

3.4. The Dirichlet problem with Holder data. We start with a lemma in the spirit of 
Theorem 3.4 in [1]. 

Lemma 3. Let D C 1R 2 be a star-like Lipschitz domain. There exist e = e(D) > 0, so that 
for 2/3 — e < p < 1 the maps 



2 
1 



+ K* : H p (dD) -> H p (dD) 
- K* : H p (dD) -> H p (dD) 



are invertible. 



Assuming the validity of Lemma 3, we can easily show part three of Theorem 2. Indeed, 
observe that -I+K : C a (dD) -> C a {dD) is the adjoint map to -I+K* : H p (dD) 



and — -I 
where a = 



- K : C a (dD) - 
1/p — 1. Thus 



C a (dD) is the adjoint map to -\l + K* : H p (dD) 



> H p (dD) 
H p (dD), 



1 

2 
1 

2 

are invertible operators for a < 



+ K : C a {dD) -> C a (dD) 
- K : C a {dD) -> C a (dD) 



1 



2/3 -e 
D a ) 



1/2 + 0(e). Thus, the solution to 



Au = 
u\ dD = fe C a {dD) 

is in C a (D). 

Remark There exists more direct arguments towards proving the C a estimates of Theorem 
2 similar to the one employed for the system of elastostatic ([2]), and later on for biharmonic 
functions ([10]). 
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The tools provided by Lemma 3 however allow for unified treatment of the problem at 
hand. More specifically, one reduces the question for solvability of the regularity and Neu- 
mann problems in H p 2/3 — e<p<2 + eto the invertibility of unitary perturbations of 
the adjoint of the double-layer potential in the same spaces. 

Lemma 3 follows from the existence and uniqueness statements for R p and N p , 2/3 — e < 
p < 1 combined with the usual duality argument. The proof of Lemma 3 is essentially 
contained in [1] (cf. Proposition 3.1-3.5). One can easily adapt the argument there to the 
two dimensional case and the extended range of p's, thus we omit the details. 

4. Existence and Uniqueness for biharmonic functions 

In this section, we briefly sketch the proofs of Theorems 5, 6. We follow closely the ad-hoc 
approach of [10], which originated in [2]. As we have mentioned earlier, a more systematic 
way of studying the problem would be the method of Lemma 3, i.e. one could build an 
operator T, whose invertibility is equivalent to the solvability of BR P and by duality to 
the Dirichlet problem with a-H61der data in the sense of Theorem 5. This program was 
implicitely carried out in [9]. We choose however the direct method for sake of clarity of the 
exposition. 

4.1. Uniqueness for biharmonic functions. An easy adaptation of Lemma 2 gives the 
following. 

Lemma 4. Suppose that for a given biharmonic function u, there is X* e D C W l , such 
that \Wu(X*)\ = 0. For p <n — 1 and p* = (n — l)p/(n — 1 — p) there is 

(15) \\M{Vu)\\ LP , (dD) < C\\M{VVu)\\ LP(dD)) 

where the constant C depends only on the distance of X* to the boundary, p and the Lipschitz 
character of dD. 

Indeed, one uses the equivalence of the area integral and the non-tangential maximal 
function for biharmonic functions as in the proof for the harmonic case to show (15). Since 
the maximum principle of [10] is valid for dimensions two and three (but not for dimensions 
bigger than three), we argue as follows. Assume that a biharmonic function solves BR P for 
2/3 — e < p < 1, such that D 2 u\ 9D = and V Tl Diu\ 9D = and \\M (Wu)\\ LP( ^ dD ^ < oo. 
Thus, after an eventual correction with a linear term, we may assume that Vu\qd = and 
u\dD = 0. By Lemma 4, we conclude 

\\M(Vu)\\ LP * {9D) < \\M(VVu)\\ LP(dD) < oo, 

where 2 — e < p* < oo. By the maximum principle, we have uniqueness for the Dirichlet 
problem in LP , hence |Vw| = 0. 

4.2. Existence for the biharmonic regularuty problem. For the existence part, we 
will use the following Cacciopolli type inequality. 

Lemma 5. Let D C W l be a domain above Lipschitz graph. Let fl± C C D be bounded 
Lipschitz domains. Let A 2 u = in D with M(V 2 u) e L 2 (dD). Let e is a small number 
as in Proposition 1. Let also 2 — e < p < 2 + e and d = dist(f2i, D \ Q 2 )- Then there is a 
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constant C , depending only on the Lipschitz constant and p, so that 



J \V 2 u\ 2 dX < C(\\Vu\\ Lp/(dDnm jM(V 2 u 



LP (3D) 



+ 



+ 



u 



LP {3D) 



+ 



+ d ^uWmnJ^Wmn^ + d 'IMIl^JI V^ILW' 

Lemma 5 appears as Lemma 5.6 in [9] for dimension three. The proof though can be easily 
adapted to this generality. We state now our main estimate for "atomic" solutions. 

Lemma 6. Let a be an unit atom in H p (dD), D C 1R 2 , with J a(z)z = 0. There exists 
e = e(D) > 0, such that the unique solution to the regularity problem 



BR n 



satisfies 
(16) 



A 2 u 

D 2 u\dD 

(V Tl ,VD lU }\ dD 
\\M(VVu)\\ LP{dD) 

J M(V 2 uf < 1, 

3D 



= 
= 

= a 

< oo 



for 2/3 - e < p < I. 

Remark The extra cancellation condition J a(z)zdz = is technical and it is possible 
to remove. Simple translation and dilation argument yields (16) for arbitrary atoms a in 
H p (dD) with the special cancelation property f a(z)zdz = 0. Since such atoms suffice to 
span H p (dD), we get 



/ 

3D 



M{v 2 u a y 



< 



\HP(3D) 



In particular, we get (16) for atoms without the extra cancellation J a(z)zdz = 0. 

Proof. (Lemma 6) We make the standard assumption that the boundary is smooth, so that 
smooth solutions exist according to the classical theory. As usual, our estimates will not 
involve the smoothness constants and after one proves the result in that fashion, a standard 
approximation technique yields (16) for general Lipschitz domains. 

d 2 

Next, observe that the boundary conditions for BR P imply that -——u(x, (fix)) = a(x) and 

dx 2 

therefore 



«<*.**)> = //««*"*»■ 



By support considerations and since J a = 0, J a(z)zdz = 0, we get u(x, tp{x)) = 0, for x > 2 
and u(x, (p(x)) = 0, for x < —2. Thus, u also satisfies a Dirichlet type boundary conditions 

A 2 u = 

D 2 u\ dD = 

x y 

u(x,<f(x)) = f f a(z)dzdy. 




Miy 2 uf I < 1. 
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The advantage of casting u as a solution to both Dirichlet and regularity type problems will 
be seen later on in the proof. 

We first estimate (16) for x-small. By Holder and L 2 regularity 

/ \ l,P / \ 1/2 

J M{V 2 uf 

\@Dn{(x,<p(x)):\x\<100} J 

For every point X G 3D fix a right cone T(X) opening upward with axis along the "time" 
axis and sides having slopes 100 || OCJ . Define the auxilliary maximal functions 

Mx(V 2 m)(X) = sup |V 2 w|, 
r(x)nr 

M 2 (VV)(X) = sup |V 2 w|. 

r(x)\r 

Mi measures the behavior of the solution away from the boundary and is somewhat easier 
to control, while M 2 captures the behavior of the solution close to the boundary. It is clear 
that M < Mi + M 2 . 

Since D 2 u = and Vw|aD £ L 2 ~ £ (dD), we deduce from the L 2 Dirichlet theory 

\\M(Vu)\\ L2 . e{dD) <l. 

Thus, interior estimates imply 

|Vw(X)| <dist(X,9D)- 1 /(2--) ? 
\VVu(X)\ < dist(X,9 J D)- 1 - 1 /( 2 - £ ). 
Therefore for \Q\ > 100, we infer 

(17) M 1 (V 2 u)(Q) < IQ^ 1 ' 1 /^ g L 2/3+0(£) (<9£). 

For R > 10 and 1 < r < 2 define the Carleson region above Zr = {(x,ip(x)) : \x\ ~ R} 

Q T = nf = {( x , t) : R/t < \x\ < Rt, ip(x) < t < IOOtH^'H^} . 
By Holder's inequality, we have 

, P/2 

72„ \2 



J Miy 2 uf < R l ' p/2 [ J Miy 2 uf 

dDndfli VDndfli 



Hence to show (16) it will be enough to prove 

M{V 2 uf < R 



ivi y v U) ^ n 2 £ , 
dDndQi 



for some positive e > 0. 

From the L 2 regularity result in Q^, we have 



J M 2 (VV«) 2 < J |VV«| 2 + J |V Ti Vm| : 

dfl T ndD dn T \dD dn r ndD 
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Averaging in r and VTiVu|an t nau — yield 

(18) J M 2 (VVu) 2 < R- 1 J \VVu\ 2 . 

dfliCidD Q 2 

The boundary conditions u\qq 2 ^qd = 0, V-u|an 2n< 9Z) = and Lemma 5 imply 

n 2 

for some eventually bigger box Q 3 C D still having diameter ~ R. Since u|dn 3 naD = 0, one 
estimates 

(20) IMIlw ~ ^11^11^(03)- 

Combining (18), (19), (20) with the obvious || V 2 u\\ L2m < R 1/2 \\M 2 (V 2 u)\\ L2{dn3ndD) yield 



(2i) /|v 2 «| 2 <it:- 1/2 ||v«|| L2(n3) J M 2 (V 2 u) 



1/2 

2„.\2 



n 2 \}n 3 ndD 



1/2 

2„.\2 



(22) J M 2 (V 2 u) 2 <R-y 2 \\Vu\\ L2m ( J M 2 {V 2 u) 
d^ndD \>n 3 ndD 

As usual, one has ||Vu|| L2(n3) < R 1/2 \\M(Vu)\\ L2{dn3ndD) < R 1/2 and by iterating (22), one 
gets for every 5 > 

(23) J M 2 (V 2 u) 2 < C s R- 2+S , 

d^ndD 

which barely fails to make f M 2 (V 2 u) p convergent. One also has 

dD 

(24) l|V 2 «L (ni) < ^||M 2( V 2 .)|| L2(ni) < CtR-Wft 
Thus in order to get a better estimate we resort to (7). Write 

r / f \ V(2-e) / r \ V(2-e)' 

(25) / |V«| 2 < / |V«M / |V M |( 2 -)' 

< i?v^v(2- £ y|| M(Vw) || L2 _ e(aD) || | VM | 4 j by (7) 

where f2 3 C Q,^ is still a domain with diameter ~ R. Since one can clearly derive (23), (24) 
for f2 4 instead of f2i (with eventually bigger constants), we use those with 5 = e/100 to 
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estimate the right hand side of (25). We get 
(26) 
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/ |V«| 2 < i? 3 / 4 -^||V.||i/ 2 ( , 4) + ^-^UVull 



L 2 (n 4 )- 



n 3 



Iterate (26) to get ||Vu|| i2 (n 3 ) < R 1 ^ 2 Going back to (22) and using the newly obtained 
bound for || Vw|| i2( -^ 3 ), we have 

1/2 

(27) 



J M 2 (V 2 u) 2 < R- 1 - ^ ( I M 2 (V 2 u) 2 
dfiindD \>n 3 ndD / 



Iterating (27) now gives the desired estimate 

(28) J M 2 (V 2 u) 2 < R- 2 -° {£ \ 

d^ndD 

which completes the proof of the Lemma. 

We finish the proof of Theorem 6, based on Lemma 6. 



□ 



Proof. (Theorem 6) Start with data (/, g) G H[{dD) x H p (dD) in BR P . Select a harmonic 
function h with Dirichlet data /. Define H = h-i to be the primitive of h. Stein's lemma 
and the regularity estimates in Theorem 2 imply 

LP(dD) < \\M(Vh)\\ 



(29) 
By Theorem 4 



Hf(dD)- 



dh 



dN 



< 



HP (3D) 



Hf(dD)- 



dh 



Lemma 6 then provides a biharmonic function v, with data g — such that 



' ^\lp{8D) ~ 



9 



dh 



dN 



(30) ||M(V 2 " 

HP(dD) 

Set u = H + v to get a biharmonic function satisfying the desired boundary conditions. 
Combining (29), (30) yields the estimate 



||M(V< 



u 



LP (3D) ~ 



(dD)xHP(dD) ' 



□ 



4.3. Holder continuity of solutions. We will show Theorem 5 based on the estimates for 
solutions of the regularity problem, in particular for the Green's function. 

Proof. (Theorem 5) Recall the representation formula for biharmonic solutions (2) 
(31) u{X) = I f (Q)^A Q G(X,Q)da(Q) + f £ f j (Q)N j (Q)A Q G(X : Q)da(Q). 
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By rescaling and dilation, it will suffice to show 



\u 



(X) - u{X*) - (X - X\ Vu{X*))\ < \\fi\\c°(8Dy 

j=i 



when dist(X,dD) = 1. To achieve that one obviously needs estimates for AG\od- 

Observe that by (28), (11) and the construction of the solution for the regularity problem 
(see the proof of Theorem 6), we can establish the following estimate for the Green's function 
of A 2 in D 



(32) 



J M(V 2 G(X,Q))dQ<2- 2k 2- 



■0{s)k 



dDD\Q~X*\~2 k 



when dist(X,dD) = 1 (cf. (3.5) in [10]). 

For the first term in (31), the normal derivatives of the harmonic function AqG(X,-) 
are converted into tangential derivatives for f and one eventually replaces the harmonic 
function AqG(X,-) by its conjugate harmonic functions (cf. [10], p. 395). The estimates 
that one needs are then the same ones as for the second term. We estimate the second term 
in (31) below. 

Note that by considering u(X) - u(X*) - (X - X*,Vu(X*)) instead of u(X), we may 
assume without loss of generality that fo(X*) = and fj(X*) = 0. Thus 



J2f j (Q)N j (Q)A Q G(x,Q)daQ 



f)D 



Y^MQ) - f J (X*))N J (Q)A Q G(X,Q)daQ 

j 

$ T,Wfj\\c« {dD) J \Q-XT\A Q G(X,Q)\da(Q) 



< 



dD 



$ Ys^c^D)^ 2 ^ ( j \^ Q G(X,Q)\ 2 da(Q) 

j k=0 \Q-X*\~2 k J 

< Vllf-ll ^ 2 {a ~ l/2 ' 0(£))k < V II Ml 

j fc>0 j 

provided a < 1/2 + 0(e). 



1/2 



<by (32) 



□ 



5. Counterexamples 

In this section, we provide counterexamples to show that the statement of Theorem 2 
(and consequently Theorem 3) is sharp. Since Theorems 5 and 6 provide an extension and 
essentially contain the results of Theorem 2, the harmonic functions considered below should 
be viewed also as biharmonic counterexamples showing the sharpness of the statements of 
Theorem 6 and Theorem 5. 
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Consider the domain 

Q s = {z e C : \z\ < 1, < argz < 2tt/(1 + 5)} , 

for some small 5 > 0. Observe that the domain has Lipschitz constant 0(1/5) and is very 
"non-convex" as 5 — > 0. We will show that the Dirichlet problem with Holder data in 
C a (dQs)i ol > (1 + 5)/ 2 is not uniquely solvable. 
Define the harmonic function 

u{z) = \m z (1+5 ^ 2 , 

with the obvious identification of the complex plane C with M 2 . It is easy to check that d£l$ 
consist of two segments and an arc, so that u vanishes on the segments and u e C l on the 
arc. Altogether, we have that u solves a Dirichlet problem on £7,5 with C 1 data, while one 
obviously cannot control more than |Mlc(i+<s)/2(n i5 )- Thus, we have shown that part three of 
Theorem 2 is sharp. 

Next, we invoke the equivalence results of Section 3 to conclude that since one cannot 
solve the Dirichlet problem with Holder data, then the regularity and Neumann problem 
must fail to be solvable as well. 
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